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A model glass is considered with one type of fast (/3-type) of processes, and one type of slow 
processes (a- type). On time-scales where the fast ones are in equilibrium, the slow ones have a 
dynamics that resembles the one of facilitated spin models. The main features are the occurrence 
of a Kauzmann transition, a Vogel-Fulcher-Tammann-Hesse behaviour for the relaxation time, an 
Adam-Gibbs relation between relaxation time and configurational entropy, and an aging regime. 
The model is such that its statics is simple and its (Monte-Carlo-type) dynamics is exactly solvable. 
The dynamics has been studied both on the approach to the Kauzmann transition and below it. 
In certain parameter regimes it is so slow that it sets out a quasi-equilibrium at a time dependent 
effective temperature. Correlation and Response functions are also computed, as well as the out of 
equilibrium Fluctuation-Dissipation Relation, showing the uniqueness of the effective temperature, 
thus giving support to the rephrasing of the problem within the framework of out of equilibrium 
thermodynamics. 

1. INTRODUCTION 

A glass can be viewed as a liquid in which a huge slowing down of the diffusive motion of the particles has destroyed 
its ability to flow on experimental time-scales. The slowing down can be expressed through the relaxation time, i. 
e. the characteristic time needed to have one inter-particle diffusion process of a particle while it is rattling between 
its neighbour particles, that form a cage around it. This relaxation time is proportional to viscosity. Cooling down 
from the liquid phase, at some point the system falls out of equilibrium: the slow liquid degrees of freedom are no 
more accessible and the relaxation time and the viscosity of the under-cooled melt grow suddenly by several orders 
of magnitude. The temperature at which this happens is defined as the glass transition temperature T g [[l]. At 
T g the heat capacity decreases in a clear way going from liquid to glassy phase and also on reheating an abrupt, 
but different change shows up. (Some universal behavior in the cooling-heating process was pointed out by one of 
us |§J|]-) Moreover, discontinuities of this kind occur also in the compressibility and the thermal expansivity. This 
looks similar to a continuous phase transition, even though the analogy is not perfect, because of the smeared nature 
of the discontinuities and because the smaller specific heat value occurs below the glass transition, rather than above, 
as would normally occur in mean field phase transitions. 

The above described transition is not a true thermodynamic phase transition, but it is strictly kinetic in origin: it 
takes place when the relaxation time becomes longer than the observation time and marks the transition from ergodic 
to non-ergodic behaviour. In general the location of this transition, the empirical glass transition temperature T g , 
depends on the cooling rate, more precisely, on the cooling scheme. The absolute glass transition temperature is 
defined as the temperature where the viscosity equals 10 13 Poise, and the equilibration time of the order of days. It is 
related to the slowest possible experiments one can realistically do. Cooling at higher rates there is a glass transition 
at a somewhat larger dynamical glass transition temperature. 

The very slow relaxation of non-ergodic systems evolving towards equilibrium structures on time scales longer than 
the characteristic time scales of the experiments, depending on the history of the system (e. g. on the phase space 
region in which the initial conditions are chosen and on the cooling rate) is the so called regime of aging dynamics 

Experimental data for the viscosity pattern of glass forming liquids are often fit to a Vogel-Fulcher-Tammann-Hesse 
(VFTH) behavior r cq ~ exp[A 7 /(T — To) 7 ], where the fitting parameter To depends on the material and the 

range of temperatures in which the fit is performed. The exponent 7 is usually set equal to 1, and an argument for this 
choice was given by Adam and Gibbs H . However also exponents 7 7^ 1 are compatible with data, merely affecting 
the width of the fitting interval. An analytic approach by Parisi |Io[ , on the other hand, using replica trick and field 
theory, gives 7 = 2 in three dimensions. Here we shall consider 7 as a model parameter, that can be chosen below, 
equal to or above unity, and investigate the aspects of this standard picture. 

Kauzmann (ll]] pointed at the paradox that the difference between the liquid entropy and the crystal entropy 
(i.e. the entropy of the most organized state for the system) if naively extrapolated to zero temperature would 
become negative at some point. To circumvent this unphysical result he proposed the occurrence of a thermodynamic 
phase transition at the temperature (commonly denoted by Tk) where this entropy difference vanishes. Such a 
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thermodynamic transition would be characterized by a discontinuity of the specific heat and by the exponential 
divergence of the relaxation time (VFTH for fragile glasses or Arrhenius for strong ones) . Connected to this last feature 
is the usual assumption that the fitting parameter To of the VFTH law coincides with the Kauzmann temperature 
Tk ■ We note that at this phase transition, the divergence of the relaxation time is not algebraic in temperature, as 
happens in ordinary continuous phase transitions, but it is exponential and no susceptibility diverges at the critical 
point. 

The residual entropy given by the difference of the entropy of the undercooled liquid and the entropy of the 
vibrational modes of the crystal that could in principle be formed, is usually called complexity or configurational 
entropy. According to standard knowledge, the Kauzmann transition should be characterized by a vanishing, or 
minimal, configurational entropy. This prediction is very difficult to test experimentally, since the relaxation time 
is too long. The existence of a Kauzmann transition was nevertheless recently supported both by analytical and 
numerical results |IJ-[l4|]. The configurational entropy is the entropy determined by the number of states that the 
system at temperature Tk < T < T g can visit. 

At a given dynamic critical temperature To, generally greater than T g , the separation of the time scales of slow 
(a) and fast (/?) processes starts to increase more rapidly than at higher temperature. Referring to the phase space 
we can say that structures get organized at two levels: some minima of the free energy are separated by very small 
barriers and between them (3 processes take place; groups of those minima are contained in bigger basins separated 
by barriers requiring a greater free energy variation to be crossed. To make the system go from the configuration in 
one of these basins to another configuration in another basin, i. e. to have an a process, a longer time is needed. The 
time scale on which these processes are happening are, however, at To and below (but above T g ), still very short in 
comparison with the observation time. In a cooling experiment the system is thus still in thermodynamic equilibrium. 
Going on with cooling, the deepness of the local and global minima, appearing in the thermodynamic potential and 
corresponding to different metastable and stable states, grows: barriers between them become higher and higher until, 
at the glass transition temperature T g , some states become impossible to reach during the time-scale we set for our 
system, i. e. the experimental time. The configurational entropy is the observable that counts the relevant states. 
As temperature decreases further the configurational entropy starts to decrease because the states available for the 
system are less and less. The Kauzmann temperature is reached when the system is stuck in one state and cannot 
move to any other, because, even asymptotically and even for short range systems where activated processes were 
present for Tk < T < T g , the free energy barriers become infinite. 

The configurational entropy density s c is usually connected directly to the relaxation time through the Adam-Gibbs 
relation ||]: r eq ~ exp (l/s c ). 

To recapitulate we have been referring basically to the following different regimes for a glass former. 

• For T > Tjj the system is in a disordered phase. Diffusion processes have a very short relaxation time. At very 
high temperature the free energy describing the system has only one global equilibrium minimum and cooling 
towards the To temperature small local minima show up. 

• Around To a dynamic transition takes place. The phase is still disordered but the free energy increases its 
roughness and some local minima become deeper: a-(3 bifurcation is qualitatively enhanced. In a simple mode 
coupling theory |l5| ] this is the temperature T mc at which a static transition is predicted with an algebraically 
diverging relaxation time. In the p-spin spin glass model Jll|[l7j this corresponds to the dynamic critical 
temperature To at which the system goes to metastable states of energy higher than the minimum energy. For 
T g < T < To the dynamics of a processes has a huge slowing down but the temperature is high enough to reach 
equilibrium on the experimental time scales. 

• Around T g , that depends on the cooling rate, another transition takes place. Many other local minima appear 
and the free energy roughness is such that the deepest local minima, corresponding to metastable states, become 
ergodically separated on the time scales of the experiment. For Tk <T < T g the system has a very slow aging 
dynamics between the metastable states, proceeding by activated processes. 

• At T — Tk a thermodynamic phase transition shows up, with exponentially diverging relaxation time. The free 
energy barriers between deep local minima increase to infinite and the system gets stuck in one single minimum 
forever. Ergodicity is broken at any time scale. In the p-spin model this corresponds to the temperature at 
which the replica symmetry is broken |]l8| . The Kauzmann temperature Tk is usually assumed to coincide with 
the fitting parameter To of the VFTH law. Below the Kauzmann temperature the system evolves only through 
the configurations belonging to the ergodic component of the phase space where the dynamics brought it during 
the cooling. 

• In a cooling experiment that goes below Tk, and in a quenching experiments to a temperature below Tk, the 
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aging dynamics visits only states with a free energy that is higher than the one in the static limit. Dynamics 
behaves as if occurring at a higher temperature. 

The exponential divergence of time scale in glasses (opposed to the algebraic divergence in standard continuous 
phase transitions) might induce an asymptotic decoupling of the time-decades. The reasonable assumption can be 
made that, in a glassy system that has aged a long time t, all processes with equilibration time much less than t are 
in equilibrium (the [3 processes), while those evolving on time scales much larger than t are still quenched, leaving 
the processes with time scale of order t (i. e. the a processes), as the only interesting ones. Indeed this assumption 
was already tested successfully in models similar to, but even more idealized than, the one we are going to discuss in 
the present paper. Those models showed a glassy regime with an Arrhenius law (rather than VFTH), like a harmonic 
oscillator model and a spherical spin model . The asymptotic decoupling of time scales that is the input for 

the present set of models could be the basic ground for a generalization of equilibrium thermodynamics to systems out 
of equilibrium || . That approach involves systems where one extra variable is needed to describe the non-equilibrium 
physics, namely the effective temperature. One of our aims will be to test this picture in an exactly solvable model 
glass; we shall see that there are domains where it does apply (namely when the VFTH exponent 7 exceeds unity) and 
where it does not apply (namely when 7 < 1). In this last case two extra variables will be needed making compulsory 
the introduction of an effective field besides the effective temperature. 

In the present paper we are going to investigate an exactly solvable model glass that shows all of the features 
that we recalled above for the much more complicated real glasses. The model is introduced in section |[ It is 
built by processes evolving on two different, well separated time scales, representing the a and [3 processes taking 
place in real glassy materials. In section ^| we introduce the dynamics that we apply to the model and we show the 
dynamic behaviour in the aging regime. We can implement the dynamics even below the Kauzmann temperature, 
thus getting insight in a regime where few analytic results are known. Even though the physics of our model is simple, 
we shall find general aspects of the results by formulating them in the thermodynamic language. The phrasing of 
the dynamic properties in terms of a generalized out-of-equilibrium thermodynamic frame is carried out in section |^, 
where we introduce effective parameters to take into account the history of the system. In section |^ we study two-time 
observables, such as correlation functions and response functions, and we look at the fluctuation dissipation ratio out 
of equilibrium |2(| . 



2. MODEL 

The model we study, that was firstly introduced in [pl|, is described by the following local Hamiltonian 

1 N N N N 

n[{xi},{Si}] = -Kj24- H J2 Xi - J J2 XiSi - L J2 Si t 2 - 1 ) 
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where N is the size of the system and {xf} and {Sf} are continuous variables, the last satisfying a spherical constraint: 
J2i Sf — N. We will call them from now on harmonic oscillators and spherical spins, respectively. K is the Hookc 
elastic constant, H is an external field acting on the harmonic oscillators, J is the coupling constant between {xf} 
and {Si} and L is the external field acting on the spherical spins. As we will see in this paper the simple local form 
of ( |2.1| ) allows us to introduce an analitycally solvable dynamics with glassy behaviour. 

In our simple model we introduce by hand a separation of time scales where the spins represent the fast modes 
and the harmonic oscillators the slow ones. Separation of time scales is one of the most important and most general 
characteristics that glasses are supposed to hold. Indeed, we assume that the {Si} evolve with time on a much shorter 
time scale than that of the harmonic oscillators. From the point of view of the spins the {xf} are quenched random 
variables and the combination Jx,i can be seen as a random field exerted on spin i. On the other hand, from the 
point of view of the motion of the {xf} the spins are just a noise. To describe the long time regime of the {xf} system 
we can average over this noise by performing the computation of the {Sf} partition function, yielding an effective 
Hamiltonian depending only on the {xi}, that will determine the dynamics of these variables. 

Summing out fast variables is a standard technique in physics. For instance, in any Landau- Ginzburg- Wilson theory 
there occur coefficients, of which the temperature dependence arises from summing out fast processes. We now do 
the same in our model. 

We perform the spin integration in the partition function using the saddle point approximation for large N and we 
get: 



Zs{{xi}) = / \ T\dSi\ exp{-0W [{*<},{#}]} 6\Y&-N\ (2.2) 
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We can define the effective Hamiltonian Ti c s{{xi]) = — Tlog Zs({xi}), obtaining 



K TN w + 

H eS ({xi}) = —m 2 N - H mi N - wN + — log — 



(2.3) 



(2.4) 



(2.5) 



This can also be written in terms of the internal energy U({xi}) and of the entropy of the equilibrium processes (i.e. 
the spins) S ep ({xi}): 



H ef f({xi}) = U({ Xi }) - TS ep ({ Xi }) 



U{{xi}) = N 
N 

S cp ({xi}) = — 



m 2 



log 



Hmi — r i 

w + T/2 
T 



(2.6) 
(2.7) 

(2.8) 



and it can indeed be verified that U is the Hamiltonian averaged over the spins and that S ep is the entropy of the 
spins. 

Another fundamental ingredient for the model is the introduction of a constraint on the phase space to avoid the 
existence of the single global minimum {xi — 0,Vi}, thus implementing a large degeneracy of the allowable lowest 
states. The constraint is taken on the {xi}, thus concerning the long time regime. It reads: 



?7i2 — m 1 > mo 



(2.9) 



where too is a fixed but arbitrary, strictly positive constant. The now obtained model glass has no crystalline state. 
This constraint applied to the harmonic oscillators dynamics is a way to reproduce the behaviour of good glass 
formers, i. e. substances for which nucleation of the crystal phase is especially unlikely even at very slow cooling 
rates (e. g. network formers B2O3 and Si0 2 , molecule organics such as glycerol and atactic polystyrene and different 
multicomponent liquid mixtures). These are substances for which there are non-crystalline packing modes for the 
particles composing them that have intrinsically low energy, thus favouring these disordered configurations. In general 
the crystal state still exists, at lower energy, but the probability of nucleating a crystal instead of a glass is negligible. 
In specific cases (binary solutions) the glassy state can even be lower in energy than the crystalline one and be 
thermodynamically stable with respect to any crystal configuration . 

As we will explain in detail in the next section, we impose a dynamics which satisfies this constraint and couples 
the otherwise non- interacting {xi} in a dynamic way. 

To shorten the notation for later purposes we define here the modified "spring constant" K and "external field" H: 



K = K 



J 



T/2 



H = H 



JL 



■T/2 



(2.10) 



We stress that K and H are actually functions of the {xi] themselves (through m\ and m 2 that occur in w). We also 
define the constant 



D = HJ + KL 



Recalling the definitions ( 2.1C ) it is useful to note that 

HJ + KL = HJ + KL 



D. 



(2.11) 



(2.12) 
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1. Statics 



The partition function of the whole system at equilibrium is: 
Z(T)= J VxVSexp H9W({a*}, {#})]* - N 



(2.13) 



dm\dmi exp < ~j3N 



K 

— m 2 — ami — w 



|log(^)-|(l + log(m 2 -m?))]} 



The new object that appears in the exponent is the configurational entropy 

!=-{! + log(m 2 - ml)) 



(2.14) 



which will be widely cons idered in section ^. It comes from the Jacobian of the transformation of variables Dx — * 
dm\drriii expZ, (see (2.3)). We can compute the large N limit of this partition using once again the saddle point 
approximation. The saddle point equations are found minimizing the function 



N 



F(T, mi, m 2 ) 



1 (K 



-m 2 — Hmi — w 



, w + T/2 , . 2 , 

log 1 - log(m 2 - ml) 



T V 2 

Denoting the saddle point values of m\ and m 2 as mi and m 2 the equations are 

#(mi,m 2 ) 



771l = 



K(m 1 ,m 2 ) 



™2= m x + - 

a(toi,to 2 ) 



(2.15) 

(2.16) 
(2.17) 



The form of the solutions mi(T), m 2 (T) is quite complicated because each of these equations is actually a fourth 
order equation, but they can be explicitely computed. In terms of the equilibrium values rnk we find the following 
expression for the equilibrium free energy density: 



F(T,mi(T),m 2 (T))=7V 



K 



m 2 — Hm\ — w (mi, m 2 ) 



TN 



' w(m 1 ,m 2 )+T/2 2 ' 
log 1 - log(m 2 - (mi) z ) 



T 



= U(T,mi,m 2 ) T S cp (T, mi, m 2 ) 

For the Hessian of /3.F(T, mi, m 2 )/iV we find the following expressions 



T J(T, mi, ma) 



H = /3 



= 0; 



o(u)+T/2) ;! ^ '(m 2 -mf) 2 2w(w+T/2) 2 ± (m2-m 2 ) 2 

\ J 3 £ _ q-i mi 

\ 2™0+T/2) 2 J (m 2 -m 2 ) 

, 2 I 2L 2 JL 

' 2w(w + T/2) 2 ^ JL £ 
= Hessian of 07i e g (mi , m 2 ) 
The determinant of the Hessian of j3F(T 1 mi, m 2 )/N, computed at equilibrium, is 

det(H) = * 3 (1 + Q^D + P„) 

2(to 2 — ) J 

that is always positive. In the formula above we introduce the abbreviations 



ra-2+ra. 



T 7 



J , T 1 

4i«(i«+T/2) ;4 ' 2 (m 2 -m^) 2 



1 



m 2 + m^ — mi 



(m 2 — m 2 ) 



2\2 



—mi 

Hessian of I(m\ , to 2 ) 



Q= - 
P 



J 2 (HJ + KL) 



K 3 w(w + T/2Y 
J (m 2 — mf) 



2Kw (w + T/2Y 



= lim Q(mi(t),m 2 (t)) = Q(m 1 ,m 2 ), 

t — >oo 



lim P(mi(t), m 2 (t)) = P(mi,m 2 ). 

£ — >OG 



(2.18) 
(2.19) 



(2.20) 

(2.21) 
(2.22) 

(2.23) 

(2.24) 
(2.25) 
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that we will often use in the following. The inverse matrix turns out to be 



m 2 



1 + Poo 2mi - 2±P C 



J' 



H = ^ L - | . (2.26) 

1 + Qoo^ + Poo ^ 2m 1 -2fP 00 4m 2 + 2(m 2 -m 2 )+4^Poo 

The elements of this matrix are the thermodynamic average of the fluctuations of mi and around their equilibrium 
values mi, m2, times a factor N, as we can immediately chek expanding F to the second order around rn\ and mi in 



(2.13). This holds for temperature high enough where, asymptotically, the constraint (2.9) plays no role 



3. ANALYTICALLY SOLVABLE MONTE CARLO DYNAMICS WITH GLASSY ASPECTS 

We assume as the dynamics a generalization of previously introduced parallel Monte Carlo dynamics for the har- 
monic oscillators. This kind of analytic Monte Carlo approach was first introduced in p3[ , and later applied in [|l9| to 
the simpler, exactly solvable harmonic oscillator model (which is just our model after setting J = L = 0) and by one of 
us J^,|| also for a spherical spin model (which is the present model after setting H = K = and considering the {xi} 
as quenched random variables). The thus obtained dynamical model with a very simple Hamiltonian and a contrived 
dynamics has the benefit of being not only programmable on a computer, but even being solvable analytically, which 
yields a much deeper insight in its properties. Moreover, in the long-time domain the dynamics looks quite reasonable 
in regard to what one might expect of any system with a VFTH-law in its statics. 

In a Monte Carlo step a random updating of the variables is performed (xi — ► x[ = Xi + Ti/yN) where the {r^} 
have a gaussian distribution with zero mean and variance A 2 . We call x (without any subscript) the energy difference 
between the new and the old state, viz. x = H({a;-}) — H({xi}). If the energy of the new configuration is higher than 
the energy of the initial configuration (x > 0) the move is accepted with a probability W(f3x) = exp(— /3x); if the new 
energy is lower (x < 0) it is accepted always (W((3x) = 1). 

The updating is parallel and it is this particular feature giving the collective behaviour leading to exponentially 
divergent time scales in a model with no interactions between particles such us ours. A sequential updating would 
not produce any glassy effect. In this sense there is an analogy with facilitated Ising models j]24(] . and with the kinetic 

25), where the transition probabilities depend 
Dehavior in situations where ordinary Glauber 



lattice-glass model with contrived dynamics of Kob and Andersen 
on the neighboring configuration; this dynamics may induce glassy 

dynamics |26| would not. Models of these types may give valuable insights in the long time dynamics, at least, within 
a class that exhibits some long-time universality. 

In a Monte Carlo step the quantities Nm\ = J2 i Xi and iVm 2 = J2i x i are updated. Let us denote their change by 
yi and y2, respectively. Following j^] we get the distribution function of y\ and 2/2, for given values of mi and m 2 ,: 

p(yi,y 2 \m 1 ,m 2 ) = J II J^2 exp ("2A 2 ) 5 x ' i ~ ^ Xl ~ yi ) 5 ~ 5^ ~ ^1 

ex P -77X2 ^27- -2T— 



4ttAV^2 ~ml V 2A2 8A 2 (m 2 -mf 
We can express the energy difference as 

K , . 

x= — y 2 -Hy 1 , (3.2) 

upon neglecting the variations of mi and m,2 that are of order (yk/N) 2 ~ A 2 /N. 

In terms of the energy difference x and of y = yi the distribution function can be formally written as the product 
of two other gaussian distributions: 

p{yi 1 y 2 \mi 1 m 2 )dy 1 dy2= p(x\m 1 ,m 2 ) p(y\x,m 1 ,m 2 ) dx dy 

- ^= exp (-^) ^= exp (JjL^m dx dy (3.3) 
\/2¥A7; y \ 2A X J \ 2A V J y v ; 

where 



G 



x = A 2 K/2, A x ^A 2 K 2 (m 2 -m 2 1 ) + A 2 K 2 (m 1 -H/Ky, (3.4) 

mi-H/K x-x A 2 (m 2 -ml) 
y{x) = , A y = — -3. (3.5) 

TU2 — m\ + I mi — H/K j r«2 — to 2 + I toi — H / K\ 

The variance of the randomly chosen updating {r{\ of the slow variables was in previous approaches [p9P,p|,p3| a 
constant. That was enough to cause an Arrhenius relaxation of the glass. To find a VFTH-like relaxation, in the 
present model we let A 2 depend on the distance to the constraint, i.e. on the whole {x{\ configuration before the 
Monte Carlo step: 

/ g \ 7 

A 2 (t) = 8[m 2 (t) - m\(t)] — ^ (3.6) 

W 1WJ \m 2 (t) - m 2 (t) - m J V ' 

where B is a constant and 7 is an exponent larger than zero that we discussed already as being used in practice to 
make the best VFTH-type fitting of the relaxation time in experiments In our model 7 is a constant; it has 

no prescribed value since we do not make any connection with a microscopic system. In the standard VFTH-law one 
would just take 7 = 1. One of our results will be to see that there are three qualitatively different regimes: 7 > 1, 
7 = 1 and < 7 < 1, showing that the situation 7 = 1 is actually non-generic. 
We also define a quantity that we shall frequently encounter in the following, 

rWE (m 2 (t)- mi (iP-m„) ' (3J) 

The nearer the system goes to the constraint (i.e. the smaller the value of m 2 — to 2 — too), the larger the variance 
becomes, thus implying almost always a refusal of the proposed updating. In this way, in the neighborhood of the 
constraint, the dynamics is very slow and goes on through very seldom but very large moves, a thing that can be 
interpreted as activated dynamics. When the constraint is reached the T becomes infinite and the system dynamics 
is stuck forever. The system does not evolve anymore towards equilibrium but it is blocked in one single ergodic 
component of the configuration space. At large enough temperatures, the combination m 2 (t) — m 2 (t) — too will 
remain strictly positive. The highest temperature, To, at which it can vanish for t — > 00, is identified with the 
Kauzmann temperature Tk- 

The question whether detailed balance is satisfied or not is also non-trivial in our model. Indeed, it happens to be 
satisfied for this kind of dynamics only for large N . For exact detailed balance we should have 

p(x\m-L, 7TI2) exp(— fix) = p(— x\mi, m 2 ) (3-8) 

but now, when we perform the inver se m ove {x^} — > {xi}, the probability distribution is also depending on 



the {ri} through A as defined in (3J3). Thus the right hand side of the detailed balance consists of a 
p(—x\m' l7 m' 2 ; A' 2 ) 7^ p(—x\mi, m 2 ; A 2 ). Expanding this probability distribution in powers of 1/N, however, we 
get that p{— x\m[, m' 2 \ A' 2 ) = p(— x\m\, m 2 ; A 2 ) + 0(A 2 /N). Other terms of 0(A 2 /N) were already neglected in the 



approximation of x done in (3.2). So, inasmuch as the whole approach is valid only for N — * 00, detailed balance 
is also satisfied; it would slightly be violated in a finite N simulation. We work at very large and, even though 
A 2 oc r(t) grows as the system approaches equilibrium (it even diverges at the Kauzmann temperature), we perform 
first the thermodynamic limit computing the dynamics equation and only eventually the limit t — > 00. If we would 
do the opposite there would be a region around the Kauzmann temperature where the detailed balance is violated 
and the dynamics is not the one discussed here. However, this is not our aim since we are interested in the ergodicity 
breaking that takes place in systems with a large number (Avogadro like) of variables. 

In the harmonic oscillator model and in the spherical spins model studied in the dynamics was performed 

within this approach, but at fixed A. Both cases showed a relaxation time diverging at low temperature with an 
Arrhenius law, typical of strong glasses. We could also study enhanced Arrhenius law by setting mo — in the present 
model but here we want, instead, to develop a model representing a fragile glass with a Kauzmann transition at a 
finite temperature. 

The Monte Carlo equations for the dynamics of mi and m 2 can now be derived according to the lines of || . They 
read: 

toi = J dyidy 2 W((3x) y 1 p{yi,y 2 \m u m 2 ) = j dxW(f3x) y{x) p(x\mi,m 2 ) , (3.9) 
to 2 = J dyidy 2 W((3x) y 2 p{y\,y 2 \m u m 2 ) = -L J dxW({3x) [x + H y{x)) p(x\mi,m 2 ) . (3.10) 
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Before performing the study of dynamics we define two new variables fj,± and /12 depending on mi and m.2 and 
representing, respectively, the deviation from the equilibrium state and the distance from the constraint: 



H 



Ml 



K 



mi 



(3.11) 

fi 2 = TO2 - mf — mo. (3-12) 
When fix — equilibrium is obtained and the equilibrium value of mi is given by the solution of the equation 

H(mi,m,2) 



K (mi, 7712) 



mi 



(3.13) 



that is the saddle point equation ( 2.1q ), with m^ being the equilibrium values of rrifc, and it can be proven equivalent 
to a fourth order equation. 

When \i2 = the constraint is reached. This will happen if the temperature is low enough (T < T ). T Q is the 
highest temperature at which the constraint is finally reached by the system. 

Above To equilibrium will be achieved without reaching the constraint. The temperature is too high for the system 
to notice that there is a constraint at all on the configurations: lim^oo ^{t) = M2CO > 0- 

At and below To the system goes to configurations that become arbitrarily close to the constraint, and then stays 
there arbitrarily long. 

When all system parameters are fixed (aging setup) the equations of motion (3.9)-( 3~10| ) become in terms of \i\ and 

M2 



fli = —JQ / dx W((3x) x p(x\mi, ra-i) — (1 + QD) I dx W{(5x) y(x) p(x\mi, 1122) 



II2 = I dx W(l3x) xp(x\mi,rri2) + 2/ii / dx W(f3x) y(x) p(a;|mi,m 2 ) 



where we have used D and Q defined respectively in (2.11) and (2.24). 
We also shorten the expression K(m 2 — mf) by the parameter 

T e = K(m 2 - mf) , 



(3.14) 
(3.15) 

(3.16) 



possibly depending on time through mi(t) and m^t). For the moment this is just an abbreviation but in the next 
section we will show that an alternative description of the dynamics is possible where T e {t) turns out to be a mapping 
of the history of the system into an effective thermodynamic parameter. This effective temperature would be the 
temperature of a system at equilibrium visiting with the same frequency the same states that the actual - out of 
equilibrium - system at temperature T is visiting on a given time-scale during its dynamics. 

In the time regime where T ^>> x 2 /T% ~ 0(1) {^{t) <C 1), the gaussian distribution of the x can be approximated 

by 



t 1 \ exp(-T) 
p(x\mi,m2) — 7= exp 



x 
2Tl 



4T e v / r 7 r 

and the equations ( 3.14 ), ( ft.lE ) become 

3(l-2r + 2r 2 ) 



1 - 



16T„ 2 r 512T„ 4 r 2 



til = 4T 



JQK(m + fi 2 )r 1- 



r 



2(m + /Lt 2 )r 1 - 



3(l-2r + 2r 2 ) 
f 



-^(l + QD) (T-(l-3r 
& (r - (1 - 3r + 4r 2 )) 



4r 2 )) 



where r is the normalized difference between the parameter T e and the heat-bath temperature T: 

T P —T 



2T P — T 



and 



T ^cxp(-T) (i _ r) 



7lT 



(3.17) 

(3.18) 
(3.19) 

(3.20) 
(3.21) 
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T (upsilon) is the leading term of the expansion of the integral representing the acceptance rate of the Monte Carlo 
dynamics 



dxW (/3x)p(x\mi , m 2 ) 



exp(-r) 



(1-r) 



l-^(l-2r + 4r 2 )+0(/i 



2-yx 



(3.22) 



The solutions to the equations ( 3.18| ) and ( 3.19| ) depend on the relative size of \i\ and fi2, thus also on 7, as well 
as on r that has a different behaviour above To, where T e tends to T in the infinite time limit, and below, where T e 
never equals the heat-bath tempe ratur e (see section^). 

The solution to the equation ( 3. 19] ) can be easily found neglecting the second term, proportional to nf. It is 
expressed in the implicit form 



erf(?T(t)) „ exp(T(t)) Snmoj, M nn . 
2ir : 2 — — = 1 + const 

i r(t)V7 7T 



where 



crf(z) = —= 

\/7T 



e"* dt. 



To second order approximation this can be written as: 



[log(t/*o) + clog(log(t/*o))] 



l/ 7 



(3.23) 



(3.24) 



(3.25) 



The constants to an d c depend on the temperature phase as will be clarified in the following. 

Above To (3.25) is more precisely the behaviour of 5fj,2{t) = fJ,2{t) ~~p2(T). Since in this range of temperature 
T e (t) — T ~ 5fi2(t), the first order expansion of r is: 



Sfi 2 (t) 



mo + fJv(T) 



2P^D 



JK eq (m +-p 2 (T))(l + Q ao D) 



(3.26) 



where and are given in (2.25)-( pT24 ). In this case in equation (3.25) c is equal to 1/2 and the expression of to 
in terms of the parameters of the model is: 



8 7 r(0)[l + P oo /(l + Q ooJ D)] 



(3.27) 



where T(0) is the initial value of T(t). 



M- 2 (t) - ^2 e ' 




10000 10000a 



exact: T(0)=10 
exact: T(0)=1 
approx: T(0)=1 





1 le+20 le+40 le+60 le+80 le+100 

FIG. 3.1. The difference between /i2(t) and its asymptotic value JT 2 (T) is plotted for heat bath temperature T = 0.41, 
slightly above the Kauzmann temperature To = 4.0 0248. The case is plotted with K = J = 1, H = L = 0.1, mo = 5. The 
upper tw o cur ves represent the exact solution (3.23) with two different initial conditions. The lower one is the approximated 
solution (3.25). In the inset the initial behaviour is shown: clearly the approximation is valid already after a few decades of 
the dynamics. 







Below To the qualitative behaviour of ^(i) (in this case the JL 2 part is zero) is t he sa me, but T is never reached. 
This implies that r goes to some asymptotic constant r^. Concerning the solution ( |3.25 ) the only difference is in the 
values 

c =l±l . tQ = B ^ (3 . 28) 

27 ' 8mo7r«,(l - re) 

In figure 3.1 we show the exact solution, numerically computed, of equation ( 3.1 9| ) for a particular choice of the 



parameter values: K = J = 1 , H = L = 0.1, mo = 5, B = 1, 7 = 2. We can see that after a couple of decades the 
behaviour is the one given in (|3.25). 



The ratio of equations ( 3.1 8| ) and (3.19) brings the equation: 

dfii _ + QD){T + 2 - 3r + 2r 2 ) - JQT e r 
diL 2 



(3.29) 



2r(m + li 2 ) - ^(T + 2 - 3r + 2r 2 ) 

With respect to the relative weight of fj,\ and 1x2 we can identify different regimes, where the solution has different 
behaviours. 

1. T > To. The leading term of the solution is give n by t he stationary solution. We can also neglect the term of 
0(n\T) in the denominator. Using the expansion (3.26) for r we get:: 

- 1 i + P. + Q-U^ (3 . 30) 



and 



Here we have also expanded T(t) as 

r W= r- 7 r^. 

M2 



(3.31) 



(3.32) 



We are most interested in what happens next to the the Kauzmann temperature, i.e. for very big T, at long 
but not extremely long times, that means 5fi2(t) small but not vanishing. A more detailed treatment, including 
an expansion in T — Tq of /J 2 appearing in P x . and , can also be done, looking carefully up to which extent 
7*2 can be approximately neglected with respect to Sfi2{t). We can neglect ~p 2 with respect to the whole 112(f) 
at temperatures very close to the Kauzmann temperature and for times that are not extremely long, so that we 
are far from the thermalization and the dynamics has still aging behaviour. In figure (3.2) we show the relative 
weight of /Z 2 on the whole Hiit) for a specific case. As is clear from the figure, as soon as we go too far from To, 
we cannot neglect with respect to ^(t) its asymptotic value Jl 2 . 




100000 le+10 le+15 le+20 le+25 le+30 

FIG. 3.2. Ratio of ~p 2 (T)/fi2 (t) at different temperatures, at and above the Kauzmann temperature. Too far away from To 
the contribution of ~p 2 (T) to fx^it) becomes relevant. The case is plotted with K — J — 1, H — L — 0.1, mo = 5. For this set 
of parameters the Kauzmann temperature turns out to be To = 4.00248 
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2. T < To, 7 > 1. In this and in the following cases the asymptotic value of (i 2 (t) is /Z 2 = so that 8n 2 (t) = 
fj- 2 {t)- Also in th is dy namic regime the adiabatic approximation can be carried out and the second term in the 
denominator of (3.29) is again negligible. In this case the leading term of r in its expansion in powers of fx 2l r^., 
is of 0(1). Therefore we get 

where 

T e = lim K(m 1 [t),m 2 (t))(m + fJ, 2 [t)) = K^mo , (3.34) 

t — >oo 

K^=K (Mi , m 2 ) , (3.35) 
T — T 

(3.36) 



2T e — T 

3. T < Tq, 7 = 1. In this case the adiabatic expansion is no more consistent. We have to solve the equation ( 3.29] ) 
taking dfj,i/dfj, 2 into account. To leading order the equation takes the form: 

dMi = MiF(l + Q D)-JQ^T e r x + + + + c(r }< (g 3?) 

Defining the quantity e = ^^^"^ we identify other five sub- regimes in the case 7 = 1. 

(a) e > 1. The solution is 

(i) = JOoTer _ ^ 1 

Zr^mo^e — lj e — 1 

The exponent e is always positive, at least in cooling, because T e > T making and positive, ci is 
also positive because it is the exponential of the integration constant (the value of which depends on the 
initial conditions). Since e > 1, the second term in the right hand side can be neglected and fi± ~ fi 2 . 

(b) e = 1. We find: 

JTer^Qn, log fi 2 (t) ( . 

= - TTqZd r(t) +C2/i2 ' (3 ' 39) 

where c 2 is the integration constant and can take any value. In the long time dynamics the logarithm term 
will take over and, independently from the initial conditions, will be [i\ > fi 2 and positive. 

(c) 1/2 < e < 1. The second term in ( 3.38| ) is leading and the solution is 

Hi{t)= Cl ^—if 2 {t). (3.40) 

ci is a positive constant and [i\ ^> fi 2 and positive. 

(d) e = 1/2. When If e < 1/2 the second term in the denominator, always neglected up to now, has to be 
taken into account. In this case the leading term in the denominator goes to zero and JTgT^Q^ can be 
neglected with respect to /J,iT(l + Q^D) in the numerator. We can thus easily solve the equation: 

dfi 2 _ 2r oo m - 2,u?r 



For e = 1/2 we get 



dfii MiT(l + Q^D) 



/fc(t) = - 1 + Q D £{t)\ogni(t)+c 2 nl(t) (3.42) 



that is not invertible analytically. It is clear anyway that in this sub-regime fi\ 3> [i 2 . c 2 can take any 
value. 
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(e) e < 1/2. The solution is 



where ci > 0. It is still \i\ 3> /J2- 
4. T < To, 7 < 1. Considering also the term /ufr in the denominator of (3.29), the solution is now: 



' ./W( i -££> )rw >- 



In this low temperature regime is \x\ 3> Hi once again. 



For 7 = 1, e < 1/2 and for 7 < 1 the solution to equation (3.41) involves only the absolute value of fii, thus giving 
two possible choices for the sign of the function ^,1(^2). In order to guarantee continuity of fi\ at the parameters 
values at which the dynam ic s cha nges r egim e, we imposed to the \x\ in two contiguous regimes to have the same sign. 



That means that in (3.42), (3.43) and (3.44) we chose the plus sign. 

The time dependent variables Hi(t), fi2(t) give the dynamic behaviour of every observable in the long, but not 
extremely long, time regime, i. e. in the aging regime. When the time increases further the dynamic will exponentially 
relax to equilibrium like exp(— t/r eq ). We will see what T eq is in the next section. 

4. OUT OF EQUILIBRIUM THERMODYNAMICS 

The history of a system that is far from equilibrium can be expressed by a number of effective parameters, like the 
effective temperature or other effective fields, in order to recast the out of equilibrium dynamics in a thermodynamic 
approach |||. The number of effective parameters needed to make such a translation is, in principle, equal to the 
number quantities considered. For a certain class of systems, however, there is some effective thermalization and the 
effective parameters pertaining to processes having the same time scale become asymptotically equal to each other in 
time. Examples of out of equilibrium regimes governed by a single effective temperature have been considered in || 
p8| . In computer glasses the approach applies with some success. p9| , |3"o| ] . 

Given the solution of the dynamics (thus the form of the functions rn\[f) and 7712 (i)) a quasi-static approach can be 
followed by computing the partition function Z e of all the macroscopically equivalent states (those having the same 
values for 7711,2) at the given time t. The measure on which this out of equilibrium partition function is carried on is 
not the Gibbs measure. In order to generalize the equilibrium thermodynamics we assume an effective temperature 
T e and an effe ctive field H e , and substitute the equilibrium measure by exp(— Tl s({xi}, T, H e )/T e ), where H e ff is 
introduced in ( |2.6| ) and the true external field H in it has been substituted by an effective field H e . T e and H e 
are at this step of the computation just fictitious parameters. However, as soon as we get the expression of the 
"thermodynamic" potential F e = —T e \ogZ e as a function of macroscopic variables mi 2 and effective parameters, we 
can fix T e and H e as taking those values that make the potential as small as possible. We thus have to minimize F e 
with respect to mi and 7712 to determine T e and H e and evaluate the resulting analytic expressions at m\ — rn\{t) 
and 7712 = rri2(t) given by the dynamics at the considered time t. Counting all the macroscopically equivalent states 
at the time t, at which the dynamical variables take values mi and m.2, we get: 

Z e (mi, m 2 ;T e ,i? e ) = J Vx exp ~^-H c s({xi},T, H e ) S(Nmi - ^ Xj) S(Nm 2 - ^ xf), (4.1) 

i i 

From this partition function we can build an effective thermodynamic potential as a function of T e and H e , besides 
of T and H, where the effective parameters depend on time through the time dependent values of mi and 7712 solutions 
of the dynamics. They actually are a way of describing the evolution in time of the system out of equilibrium. The 
effective free energy takes the form: 

F e (t) =U( mi (t), m 2 (t) ) - TS ep (mi (t) , m 2 (t) ) - T e {t)l (mi (t) , m 2 (t) ) + [H- H e {t))Nm x (t), (4.2) 

with 

T e (t) = K(t)(m + fx 2 (t)), (4.3) 

H e {t) = H -Km(t) . (4.4) 
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where the last term of F e replaces the —HNmi occurring in U (see eq. (2.7)) by —H e Nmi, and where 



N 

Z(*) = y{l 



log [m + n 2 (t)}} 



(4.5) 



is the configurational entropy and where U and S ep are given in (2.7) and ( 2.. 

As we see from (4.3) and (4.4) in the dynamic regimes 1 and 2, reported in section where fix <C (i2, the effective 
temperature alone is enough for a complete thermodynamic description of the dominant physic phenomena (H e = H) , 
while in the regimes 3a, 36 (p,\ ~ ^1%) and in 3c, 3d, 3e and 4 (/ii 3> ^2), when /ii becomes no more negligible, the 
effective field H e is also needed. 



a. Effective temperature from generalized first law 



Calling M = Nmi, and using (2/7), ( p.q ) and ( 2.14 ), the differential of the free energy (4.2) turns out to be: 

dF = -S cp dT - 2dT e - MdH e , (4.6) 

thus implying 

dU = TdS cp + T e dl + (H e - H)dM - MdH . (4.7) 

Using this expression we are able to write down the first law of thermodynamics dU — &Q+&W, in the two temperature- 
two fields case, where the change in work done on the system is, SW — —MdH. In order for the conservation of 
energy to be satisfied the heat variation has, then, to take the form 

dQ = TdS cp + T e dJ + (H e - H)dM . (4.8) 

This is the same expression obtained in the two temperature picture of |28| where the fields where absent. At equilib- 
rium, where H e — H and T e — T, this reduces to the usual expression for ideal reversible quasi-static transformations 
&Q = Td S, w ith the total entropy S = S cp + X. 



From (4.8) the complete expression for the rate of change of the heat of the system turns out to be: 



The heat flowing out of the system is — Q. Referring to the aging regimes described in section [|the quantity Q turns 
out to be proportional to fi2 in the regimes 1 (T > Tq) and 2 (T < To, 7 > 1). In the dynamic regimes 3a and 36 
(T < T , 7 = l,e > 1) is Q oc fii+fi 2 - For 3c, 3d, 3e (T < T , 7 = 1, e < 1) and for regime 4 (T < T , 7 < 1) Q oc fa. 

In every dynamic regime fa and fa are negative and this implies that the heat flow of the out of equilibrium system 
is positive in its approach to equilibrium, as it should, no matter the values of the parameters of the model. 

Starting from the first law of thermodynamics, we can derive the effective temperature in yet another way, through 
a generalization of the Maxwell relation T = dU / dS valid at equilibrium for a system of internal energy U and entropy 
S, with the derivative taken at constant magnetization (or volume). We put for simplicity H e = H in the rest of this 
subsection. Out of equilibrium, together with the previous Maxwell relation for equilibrium processes (where S has 
to be substituted by 5 cp ) also the following generalization holds: 



T _9U 
e ~ dl 



(4.10) 



A more feasible identity, where the variable to be kept constant during the transformation is the bath temperature, 
rather than the entropy of the fast processes, can be obtained pl| , p2 12 1. Let's introduce with this aim the function 

$ = T e +T e T (4.11) 

inducing d$ = T e dI—S cp dT. Through this auxiliary potential function $ we can then rewrite the effective temperature 
as 
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Te(t) 



ax 



(4.12) 



This result is a firm prediction for systems that satisfy the assumption of a two temperature the rmod ynam ics. For 
underlying mechanisms in specific cases see |3l] , |32| , |l2| . Writing the latter as <j>/Z and using ( fo.9| ) and ( 3.1C ) we get 
(neglecting terms of O(fii)): 



in agreement with (4.3) 



T e {t) = k{ mi {t),m 2 {t)) [m 2 (t) - m?(t)] 



(4.13) 



1. Statics 



To is defined as the temperature at which the constraint is reached from above: some configurations become 
infeasible and the valleys of the free energy landscape are divided by infinite barriers. The breaking of the ergodicity 
in a landscape with many minima gives rise to a real thermodynamic phase transition |TT| . 

When the constraint ( |2.9| ) on the phase space of the {x{\ is first reached, at To, in the infinite time limit, T goes to 
its minimal value T = TjT ) = 1 + logTO . Zero configurational entropy would mean that only one configuration is 
allowed for the system. Coming from high temperature there would thus be a transition from a many (meta-stable) 
states phase to a phase in which the system is stuck forever in one single minimum. This transition is what is thought 
to happen in real glasses, at the so called Kauzmann temperature. Since we are using the continuous variable {xi}, the 
entropy T (as well as S cp ) is, in our case, ill defined at low temperature: it would diverge like logT at zero temperature 
if no constraint would be present. Our value 1(Tq) is greater than zero, because this entropy counts all the mul tipl e 
ways in which the continuous harmonic oscillators can arrange themselves in order to satisfy the constraint ( |2.S| ). 
Since we are dealing with classical variables we can bypass this inconvenience just subtracting from T the constant To 
to make I(T = To) = 0. The entropy value To is related to dynamics on time scales where all the degenerate minima 
are sampled. These are much longer than the scales of our interest, and for our purposes the constant To plays no 
role. 

To see how the transition takes place we first look at the asymptotic behaviour of the effective temperature. When 
T > T and t — > oo, T e becomes the heat bath temperature T. When T < T , instead, T e neve r re aches such a 
temperature. It rather goes towards some limit value T e (T) that we can get from the equation (4.3). It may be 
rewritten for clearness in the explicit form 



(Kmo -T e + T){Kmo - T e )(JT e ) 2 + D 2 mo{Kmo - T e ) 2 - J 4 m (T e ) 2 







(4.14) 

T = Tn 



a quartic equation for the effective temperature in the infinite time limit. The same equation evaluated at T e 
gives us the value of the Kauzmann temperature To as a function of the parameter of the model. In figure (4.1) we 
plot T e versus T for a choice of parameter values. 




1 2 3 4 5 6 

FIG. 4.1. In the static regime the effective temperature is shown as a function of the heat bath temperature. At high 
temperature they coincide but below the Kauzmann transition T e does not reach T, not even in the infinite time limit. The 
system remains out of equilibrium for ever. Values of the constants are K — J — 1, H = L = 0.1, mo = 5 
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From (4.14) or Figure 4.1 we observe that dT e /dT\ T - < 1 whereas, coming from above the Kauzmann temperature, 
one has, of course, dT e /dT\ T + = 1. The derivative of T e (T) shows thus a discontinuity at T — Tq. Any thermody- 
namic function, like U and m±, will depend on the heat bath temperature both explicitly and through this effective 
temperature. For the specific heat we will have, for instance: 



C 



1 dU 




1 dU 




1 dU 




N~dT 


H 


NOT 


H 


' Ndf~ e 





(4.15) 



H 



c 2 (dT e /dT) assumed originally by Tool |33| for the study of caloric behaviour in 



This is of the same form C = c\ 
the glass formation region. 

The discontinuity in dT e 
called magnetizability in 

of these quantities contains terms proportional to dT e /dT\ H . One could now discuss the Ehrenfest relations between 
these discontinuities, and the Prigogine-Defay ratio, as was done for related models by one of us (2gJ^||. Because of 
the close analogy between all these cases, we shall not go deeper into this at the present moment. 



■,ldT\ H causes a discontinuity in the specific heat and also in the quantity — dmi/dT\ H , 
J,|l (it is the analogue of a thermal expansivity for the model here described) , because both 



2. Dynamics 

The relaxation time is the characteristic time on which the system initially out of equilibrium (because, for instance, 
of a sudden quench to low temperature) relaxes towards equilibrium. It can be defined, for instance, from the dynamical 
equations of the internal energy per harmonic oscillator u = U /N 



(4.16) 



' eq 



or, equivalently, from the equations of motion for mi, m 2 as the time at which the quantity of interest goes to 1/e of 
its initial value. In any temperature regime comes out that the relaxation time has an exponential behaviour in T: 



exp 



-Y 



(4.17) 



Making use of the solution ( 3.25 ) we find the following behaviour for the relaxation time versus the heat bath 
temperature: 

1. T > Tq. jU2(t) — ► Jl 2 {T) and near enough to the Kauzmann temperature we can linearize the latter in T — Tq. 
For t — > oo we get an exponential decay with relaxation time 



r eq oc exp 



A = B 



A 



T-T 



dT 



BK X (K - iL)(l + DQ^ + PJ) 



(K - XJ(1 + DQJ - K X P X 



(4.18) 
(4.19) 



This behaviour is a generalized Vogel-Fulcher-Tammann-Hesse (VFTH) law p-p|, where 7 can have any value 
and in particular 7=1. Looking at the configurational entropy, since at the first order expansion in fj,2 we have 
J — Jq ~ 2^~M2, w e also find from (4.17) the Adam-Gibbs relation M: 



T eq oc exp 



NB 



2mo (1 — Iq 



(4.20) 



Far from equilibrium, in the aging regime where the relaxation is very slow, we can still define a time dependent 
"relaxation time" giving the characteristic time scale on which the a processes are taking place. Always for T 
very near to Tq, in the aging regime, JL 2 , the static part of fj, 2 , is negligible with respect to the dynamic part b\i 2 
so that for the effective temperature we have the following expansion: 



T e (t) ~ T 



1 + 



-K x dfi 2 {t)+0{T-T Q ) 



(4.21) 
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we get 



where ^4(To) < j4o, meaning that in the static r is more divergent. 

2. T < T ,7 > 1. For T < T the relaxation time always diverges for t — > oo. However, as it was done in the 
case above To for the relaxation in the aging regime, an instantaneous relaxation time can be considered and 
expressed in terms of the effective temperature using the first order expansion of T e in /j,2- 

T e (t) =T e + 1 \ D ®~+ D P ~ K oa » 2 (t). (4.24) 



We find 



^-Hr.r.w)^! I (4-25) 



where A(T) is the one wr itten in (4.23). The aging behaviour just above and well below T are thus intimately 
related. The expression ( 4.25 ) resembles a VFTH law where the heat bath temperature has been substituted 
by a time dependent effective temperature T e (t) and the Kauzmann temperature by the asymptotic value T e . 
Such a relation for the time scale of the aging dynamics could hold very well in more general systems. 

3. T < To, 7 < 1. In these regimes, where /zi cannot be neglected with respect to ^2, there is no simple expression 
for r. 



5. TWO-TIME VARIABLES: BREAKING OF TIME-TRANSLATION INVARIANCE AND THE 

FLUCTUATION-DISSIPATION RELATION 

In this section we compute the correlation and response functions that, unlike the energy and the quantities mi(t) 
and m-zit), depend in a non-trivial way on two times, when the system is out of equilibrium, thus showing directly 
the loss of Time Translation Invariance (TTI) with respect to the case at equilibrium. The aim of computing such 
quantities is also to build a Fluctuation Dissipation relation and look at the meaning of the Fluctuation-Dissipation 
Ratio (FDR), d v C{t,t')/G{t,t'), far from equilibrium. 

The correlation functions between the thermodynamic fluctuation of a quantity m a (t) at time t and that of a 
quantity m b (t') at a different time t' are defined like: 

C ab (t,t') = N {6m a (t)6m b (t')} , a, 6 =1,2 (5.1) 

where (....) is the average over the dynamic processes, i. e. the harmonic oscillators. 

The response of an observable m a at time t to a perturbation in a conjugate field H b at some previous time t' takes 
the form: 

G ^'*') = ^gfp M = l,2 (5.2) 

In our model H\ = H and H2—K/2. 

Since we will very often make use of the derivatives with respect to mi and TO2 of the integrals given by the Monte 
Carlo dynamics introduced in section pi we show them in appendix explicitly computed and we shorten the notation 



defining the variables fk and gu, k — 0, 1, 2, 3 in (A.11)-(A.18) 
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1. High temperature case: T > To 



First we analyze the case above the Kauzmann temperature. In this case the expansion of the ( A.11 )-( A.18 ) in 
powers of H2{t) becomes both an expansion in S(X2(t) and in ~p 2 (or equivalently in 1/T = (JI 2 /B) 7 ), because we are 
interested in studying what happens for large times and near the Kauzmann temperature To, i.e. for small values 
of <J/i2(i)and for small values of ~p 2 (or large values of V). In the and gk written in appendix we left the notation 
without the bar meaning that an approximation for long times, i.e. an expansion in S/i2(t), or an expansion in ~p 2 , 
has still to be done, depending on the kind of approximation that we need. For the sake of clearness we repeat here 
the expansion of T(t) in the aging regime: 



The following exact relations hold: 

dmiMi = -1 - LQK, 



r(*) = r- 7 r- M 



/'2 



J 



QK 



(5.3) 



(5.4) 



We stress that <9 mi /ii and d m2 \j,\ are still functions of /i 2 , through Q and K, and that they can thus be expanded 
in powers of <5yLt 2 , leading to corrections to the fk and gk- Since in the end they will only appear in the combinations 
d mi l^i + 2m\d m2 \M\ and d mi ^i — 2L / 'Jd m2 fj,i we just give the expressions of these ones: 



d mi fii + 2m 1 d„ l2 Li 1 = -(1 + QD) 

<9miMl ~ 2— -<9 m2 /il 



(l + Q x D)-6 t i 2 Q 1 D, 



J 



1, 



where 



Qi 



Go 



J 2 {m +-p2)(3w + T/2) 



(l + Q !XI £))(mo+M2) 
is the coefficient of <5^t 2 (i) in the first order expansion of Q: 

Q{t) = Q, x + Qi6fi 2 (t) 



2w 2 (w + T/2) 



3P X 



(5.5) 
(5.6) 

(5.7) 

(5.8) 



In the following formulae the derivatives of [A\, as well as /iti itself, have to be considered as general, regular functions 

Of [12- 

Besides the expansion (5.8) of Q we also give the expansions to first order in 5fj,2(t) of the quantities r (defined in 
(pOl)), mi, K ( pl0| ), p (^25| ) and Q &2$ ): 



'(*) = 



1 



m + fi 2 
K(t) = K, 



1 + Q^D + P x 
K^P^ 



<W*) 



m\{t) = mi 
P{t) =P oa -t 



(l + Q oo L>)(m +7I 2 ) 
P~£ 

K oo J{l + Q oc D)(m +Jl2) 
P„ 



Sfjn(t) 



(1 



,D)(m + fi 2 ) 



1 + Q^D - P x 



J 2 (m +-p 2 ){3w + T/2) 



2w 2 {w + T/2) 



Sfi 2 (t) 



(5.9) 
(5.10) 
(5.11) 
(5.12) 



where w = y/ J 2 rri2 + 2JLm\ 
this dynamic regime 



L 2 + T 2 /4. For what concerns the terms containing T/ii, from ( 3.31 ) we saw that in 
JQJT 1 + Q^D + P^ 



I>i 



m +Ji 2 {1 + Q^D) 2 



-5/J.2 + O (Sfi 



(5.13) 



In [|3| equations of motion for simpler models were obtained. Our present model share the bas ic attributes that are 
needed to get those equations, namely the possibility of writing the transition probability (3.1) of the Monte Carlo 
dynamics as the product of two gaussian probability distributions (3.3), functions respectively of the energy variation 
x and of the variation y of the magnetization like quantity Xi. 
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We thus recall here the following equations holding for the equal time correlation functions: 

a+b-2 



jC ab {t,t) = J W(px) y a (x)y b (x) + A y 



E 

c=l,2 



_0_ 

dm 



[y a (x)C c b(t,t) + y b (x)C ca (t,t)} \ p(x\m ll m 2 ) dx, (5.14) 



0,6 = 1,2. 



where (recalling (|3.2|), (3.5) and using (3.6) and (|4.3|)) 

Ml 



2/i (aO 

1/2 0) 



X — £ 



77l 2 - ™? + Ml if 



41>i - mi — 



A 2 if 



x + H y r (x) 



= 4T e r 

2 

A 2 (to 2 - to 2 ) 
7712 - mf + fjif 



8(m + M2 )r + o(r M 2 ) 



(5.15) 
(5.16) 
(5.17) 
(5.18) 



Expanding the integrals J W((3x)y a (x)y b (x)p(x\mi , m 2 )dx and Aj, J W({3x)p(x\mi, m 2 )dx up to order T, the equa- 
tions become: 



C 11 (t,t) = 8[m + n 2 (t)}T(t) 1 



1 - 2r + 4r 2 



C l2 {t, t) = Wm x {t) [m + iMi(t)] r(t) ^1 - 
C 22 (t, t) = 32m 1 (t) 2 [to + /i 2 (i)] r(t) ( 1 



r(i) 

1 - 2r + 4r 2 



r(t) 

1 - 2r + 4r 2 



T(t) + 2/i(t)Cn(t, t) + 2 9l {t)C 12 (t, t) 
T(t) + 



(5.19) 



r(t) 



(5.20) 

+/a(t)Cn(t, t) + [/i(t) + <&(*)] C 12 (i, t) + 0l (t)C 22 (t, t) , 
T(t) + 32[TO + M2(i)] 2 T(t)+ (5.21) 

+2/ 2 (t)Ci 2 (t,t) + 2ff2(t)(722(t,t) • 



Due to the complicated form of the equations we are not able to find solutions valid at every time. We are obliged to 
find approximate solutions valid on given time scales. First we will study the solutions in the aging regime, for times 
that are long but not longer than some given time-scale, t g , after which the system begins to thermalize. Afterwards 
we will study the dynamics of correlation and response functions for times much longer than t g , while the system is 
approaching to equilibrium. 



a. Dynamics in the aging regime 

In the aging regime, for temperature just abo v e the Kauzmann temperature To, we can neglect ~p 2 with respect 
to Sfj, 2 (t). This means that in expressions (5.8)-(5.12) we have to put ~p, 2 equal to zero everywhere, including the 
constants Q^, P x and K x (defi n ed res pectively in (2.24), (2.21) and ( |3 .35 )) and we can write <5/it 2 (i) = fJ- 2 (t). 

To find the solutions to (5.19)-(5.21) we can firstly perform an adiabatic expansion neglecting the time derivatives 
of the correlations functions. Indeed, to first order of approximation C a b is proportional to fi 2 : they are of 0(<5/i 2 T), 
negligible with respect to the right hand side terms. We then compute the second order corrections. The solutions 
for the case T > To, in the aging regime with negligible ~p 2 and r proportional to [J> 2 (t), turn out to be: 



C n (t,t) = 
C 12 (t,t) = 
Oh (*,*) = 



1 


1 - 




D 




l 




1 - 




L D 




1 




1 - 


\- Q a 


: D 



1 - 



to + n 2 (t) 
2mi(t)m + fi 2 (t) 
4toi(£) 2 to + fj, 2 (t) 



m QiD 



2mi(i) 1 - 



m QiD 

T+qZd 



moQ^D 



M*) 2 (i-^4)--i''v»,(.)J) 



i 



(5.22) 
(5.23) 



+ 0[-)+0(p 2 2 {t))\ (5.24) 



1 + Q^D 

To get these expressions it is enough to keep in /i, g\, f 2 , g 2 (defined in ( A.13 )-( A. 16; )) only terms up to 0(T) 
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Once we have the equal time solutions we can solve the equations for the two times functions. Always following 
the approach of J^] we get the equations: 



d t C ab {t,t') = f a (t)C lb {t,t') + g a {t)C 2b (t,t') a, 6= 1,2. 

where f a and g a are defined in ( |A.13 )-( A.16| ). We introduce the function 



f = h + 



2mi5i f 2 - 2toi/i 

3i 

g 2 - 2m x g\ g 2 - 2mi#i 



(5.25) 



(5.26) 



= -4T <^ (l + QD)r- 



l + QD- 



2DQP DP(l + QD) 



l + QD j(l + P + QD) 



O 



where rhi is obtained from (3.9) as 

mi ~ 4/xiTr - ViT(l - 3r + 4r 2 ) = C(^ 2 T) 

and is negligible with respect to the leading orders. 

The decoupled equations for Cn and C\ 2 that we get are, in this notation: 



(5.27) 



(5.28) 



d t C lb (t,t') =f(t) C lb (t,t'), b= 1,2 . 

To the leading order the two correlation functions above are connected to C 2 \ and C 22 in the following way: 

C 2b {t,t')~2m 1 {t)C lb {t,t'), 6=1,2. (5.29) 

Defining the time evolution function for the considered time-scale sector as 



h(r)=exp - / f(t)dt) 



the solution of (5.28) comes out to be 



C lb (t,t') = C u (t',t') f ^- + 0(^T) 
h(t) 



(5.30) 



(5.31) 



Following the approach of [|| we also derive the response function. Neglecting the terms of C(T 2 ) (called switch 
terms in |3j) they are: 



Gn(t,t+) = -(3 / dyxdy 2 W'{l3x)y{p{yi,y 2 \m 1 ,m2)=-l3 / dxW'{j3x) [y^x) 2 + A v ] p(x\m u m 2 ) 



G 12 (t,P 



4Tr 2T 

-13 I dxW'(f3x)A y p(x\ mi ,m 2 ) + O^jT) ~ ^ =- + 0{n 2 T) 

K K 



^Tr-^ + 0( M2 T) 
K K 



G 22 (M + 



-/? J dy 1 dy 2 W'((3x)y 1 y 2 p(y 1 ,y 2 \m 1 ,m 2 ) = 

J dxW'((3x) [xy^x) +H [y^xf + A y ]} p(x\ mi , m 2 ) 
-20-^ J dxW ((3x)A y p(x\m 1 , m 2 ) + O (J^-T 
-0 J dy l dy 2 W'{f3x)ylp{y 1 ,y 2 \m 1 ,m 2 ) = 

J dxW'(px) {x 2 + 2Hxy 1 (x) + H 2 [y^xf + A v ] } p(x\m u m 2 ) 
-4/3-^ J dxW'(px) (x 2 + H 2 A y ) p( 



(5.32) 



(5.33) 



x\mi,m 2 ) + O ( =-T 



16m? 8m 2 T 



K 



-TV 



K 



(5.34) 



32Tmg + O (/i 2 T) . 



The equations describin g the evolution in t of the response to a perturbation at t' have the same shape of those for 
the correlation functions (5.25). The solutions are then 
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G a b(t,t') — G a b{t', t') S - . 

hit) 



(5.35) 



With these results we can generalize the Fluctuation Dissipation Theorem (FDT) defining another effective tem- 
perature, Tf D , by means of the ratio between the derivative with respect to the initial time (also called "waiting" 
time) t' of the correlation function C\\ and the response function G\\. 



rpFDl, ./% _ 9 t 'Cn(t,t') 

e ( ' j " Gn(M') ' 



To compute it we need: 



cVCn(M') = (^C u (f,f)) - /(*')^Cn(i',*') ~ -/(O^Cn^.f) 

h(t) h(t) h(t) 



(5.36) 



(5.37) 



-4T(0 {(1 + Q^D + Q 1 D l x 2 {t'))T{t') 

m QiD 



1 



1 + QooD 



mo + 1 - 



1 + Q„£> 



0(1)} 
hit) 



-4T(O{[m + ^(f)]r(O + O(l)} 



Eventually we get 



T e FD (<,0^ ?;(*') 



1 + 



/ i 



Vr(t') 



<?(Mi') 2 ) 



= T™(i') 



(5.38) 



where T e was first introduced in section |^ and later on derived in ([D^). We recall that T^ 1 c>c fi^. As we see here 
the above defined fluctuation dissipation effective temperature coincides, in the time-scale of our interest, with the 
effective temperature T e that we got by the quasi static app roach , only if 1/r is negligible with respect to [i 2 . This 
is true only if 7, the exponent of the generalized VFTH law ( 4.18 ), is greater than one. Otherwise the last correction 
is no more sub-leading: already for 7=1, Tf D — > T e only in the infinite time limit, i. e. for time-scales longer 
than those of the considered aging regime. As already discussed in section ||, where we presented the results of the 
dynamics of the one-time observables, the value of the exponent 7 discriminates between different regimes . For 
7 > 1 an out of equilibrium thermodynamics can be built in terms of a single additional effective parameter (the 
effective temperature T e ). For 7 < 1, T e alone does not give consistent results in the generalization of the equilibrium 
properties to the non equilibrium case and in order to cure this inconsistency more effective parameters are probably 
needed. This discrepancy was clear, from section |3|, for the regimes below To where already the one-time variables 
had different behaviours depending on the value of 7 being greater, equal to or lesser than 1. For T > Tq there was 
not such a difference at at the one-time level. It shows up, instead, at the level of two-time observables, as we just 



6. Approach to equilibrium 



For times longer than the aging regime time-s cales t he te rms that are relevant in the equations ( 5.19 )- (|5.2l| ) for 
the correlation functions and in the expressions (5.33)-( p735 ) for the response functions are different. When t 3> t g 
the equilibrium value ~p 2 of the variable /l«2 is no more negligible with respect to its time dependent part 8/12 it) (that 
i s eve n tually going to zero as t — > 00). We are in a regime where r ~ (T e ~ T). In the resolution of the equations 
(P^-(j5~2ll) this means that all the terms of OOrX^f/ 1^2) = 0(8/j,2T) are now sub-dominant with respect to those of 
0(TT) and of 0(T). To solve equations ( 5.19| )-(5.21) we use the adiabatic expansion, as in the previous case. 

The solutions for very large time, with finite, though small JL 2 and a vanishing r, are: 



1 + QD + P 
toq + /I 2 



1 



-(I + Pee) + c n< s m 8fi 2 (t) 



C 12 (t,t) = C 2l {t,t) * ff+ffij, (2rh.it) - 2§P 



Cl2,r r 



(5.39) 



(5.40) 
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c 22 (t,t) 



l + Q x D + P 
m + fi 2 (t) 
1 + QD + P 
m + 7I 2 



Cl2,5 M2 Sfj, 2 (t) 



Where 



Cll,r 



Cl2,r 



C 22l r 



1 + Q X D + P 

A 1 (m + ^ 2 (t)fK 2 J 2 Q 2 
' fi 2 (t)(l + DQ)(1 + P + DQ) 2 

^(m + ^ 2 (t)) 3 KJQ 



4mi(t) 2 + 2 (mo + M*)) + 



4to 2 + 2(to +7I 2 )+4— P c 



C22,r r 



(5.41) 



c 22 ,Sfj. 2 5jj, 2 {t). 



(1 + P) 



^ 2 (l + DQ)(l + P + DQ) 2 
16 7 (TO + /i2W) 3 KJQ 



(1 + P)(l + KLQ) + KJQ [ P- - mi 



L 



P- - mi 1 + A'LQ 



and 

Cll : d> 2 
Cl2,(5/i2 



1 



M 2 (t)(l + ^Q)(1 + P + £>Q) 2 V' J 

PiPO^ + Qi^D(l + Po) 



1 + DQ„ + P x 



1 + PQ^ + P x 



1 + P*, - (too + /Z 2 ) 



1 + P. + DQ, 



8 7 DQ OC P (1 + Po)(to + /I 2 ) 



L ( D 

mi - -P^ 1 H = 

J \ LK^l + DQ^) 



(5.42) 
(5.43) 
(5.44) 

(5.45) 
(5.46) 



(to +JI 2 ) [Pi (mi + §(1 + DQJ) + DQi (tox - ^)] 



4 7 £>P« 



1 + P x + DQ^ 
(1 + P oc )(l + K^LQJ) + K^JQ^L/JP^ - mi) 



JA^l + PQJ 2 



C22,^t 2 



1 + PQ^ + P. 



2(to + 7I 2 ) + 2to 2 + 2— P 



4toiP 00 P» 

S^j(i + dqI) 



(5.47) 



(to + 7Z 2 ) [Pi (2to 2 + 2L 2 /J 2 (P oa - 1) + too + /Z 2 ) + PQi ( 2 ™i + 2L 2 /J 2 P ao + m Q + p 2 ) 



1 + P^ + £>Q„ 



8 7 ^P C (mi - L/JPoc) (l + K ooLQ, 
JiL/^l + PQoo) 2 



where Q\ is defined in (5/7) and Pi is the first order coefficient of P of its expansion in 5fx 2 (see ( 5.12 )). 



To find the solutions to first order approximation in 5[i 2 (t) it is enough to keep in each of the ( A.13)-( A. 16) only 
terms up to 0(rT) and O (rT / (I> 2 )). Whether or not the terms of 0(rT) or those of O (rT / (r> 2 ) )7"with n = 1, 2, ... 
are the most important depends on the value of VFTH exponent 7 : if 7 > 1, O (rTfj, 2 T) <C 0(rT), for large times; 
they are of the same order at 7 = 1, while if 1/2 < 7 < 1, O (rT / (T[i 2 )) ^> 0(rT). Furthermore if 7 < 1/2 also 
terms of order O (rT / (r 2 /i 2 )) will be more important or as important as those of 0(T). For yet smaller values of 7 
more and more terms of th e kin d rT/(T n fi 2 ) will be much greater than 0(rT) in the aging regime. 

As t — > 00 the solutions ( ^3.39[ )-( p^42 ) coincides with the elements of the matrix ( 2.26 ), the inverse of the Hessian of 
the free energy of the model, i. e. they coincides with the average sq uared fluctuations at equilibrium. 
Once we have the equal time solutions we can solve the equations ( 5.25 ) for the two times functions. 
The function / is now: 



f = h + 



2rriig 1 



92 

-4T(r - 



2toi<7i 
1 + QD 



.91- 



1 



f 2 - 2m 1 f 1 
2toi<7i 

.r 2 



92 - 
p 



(5.48) 



1 + P 



4T— r 7 (m + fx 2 ) 
M2 (1 



QDP 



py 



o 



fTTr 

V M2 
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The decoupled equations for Cn(t,t') and Ci2(t,t') that we get are always the (5.28). From them we can compute 
C2i(t,t') and C 22 (t,t') as follows: 



r „ 2mi(t) 2 + 2P(t)L 2 /J 2 + mp + ^ 2 {t) 
C22 ^ t ] ~ mi (t)-P(t)L/J Cl2{t ' t] 

Using the time evolution function h for the considered time-scale sector as 

h(r) = exp f - / f(t)dt 



the solution of (|5.28|) comes out to be 



h(t) 



(5.49) 
(5.50) 

(5.51) 

(5.52) 



In the leadi ng te r ms of our expansion in d/i 2 (t) and fj, 2 the expressions for the /'s and g's are given, for the case 

1 b y dOah-dng ), 



Using ( |5.39D , fl5.40j ) and ( p.4gQ , we get 
Cu(t, * 1 + Q + ^ F No + M 2 + <9(<W))] cxp |-4^ T(t")(T(t") - l) l±gg^±p!) d t"} (5.53) 
Cu& - Y;Tv^ P \ No + M 2 + 0(^2(*))] cxp {-4 /* T(0(r(0 - l) 1 + ^ + P(O <**"} ( 5 - 54 ) 



2. Low temperature case: T < To 

Our approach allows also to study the regime below below the Kauzmann temperature To. In this last case, though, 
we have qualitatively different behaviours depend ing o n the value of 7, i.e. on the relative weight of /ii and \x 2 . We 
describe here the case 7 > 1, where fix <C fi 2 (see ( |3 .31 )). For 7 > 1, according to the results shown in section ||, it is, 
indeed, not necessary to introduce any effective thermodynamic parameter other than the effective temperature, and 
the analysis can be carried out in a way similar to the one of the previous case. In expanding the time dependent 
coefficients of C ab in the equations of motion (/x,2 and 51,2 given in equations ( |A. 13 )-( A. 16) ) we have now to take 
into account that r never vanishes, while on the contrary the asymptotic value of /i 2 {t), denoted by ~p 2 , is zero. The 
leading terms in the fk and gk (k = 0,1,2,3) are in this case those of 0(TT/^l 2 ). The sub-leading terms are those 
of order TT and T/i 2 T (coming always with r as a multiplicative factor) and those of order T. All these terms are 
diverging terms, in the limit t — > 00, and hence terms of 0(1) are now negligible with respect to them. They would 
lead to corrections to the FDR of order 1/r <C fJ, 2 ■ 



The equations of motions for the equal time correlation functions are identical to the ( 5.19 )-( 5.21 ). What change 
are the time dependent coefficients f\,2 and 51,2 and r — + 0(<5/i 2 (t)), where is defined in (3.36). 

Solutions to these equations are obtained, as before, in the adiabatic approximation and expanding all the functions 
in powers of 5fi 2 (t) and JL 2 : 



C ab {t) = [2 TOl (t)} 



a+b-2 



1 



m + fi 2 (t) 1 - 



TTqId 



mo 

r 



a\r + a 2 r 
(1 + Q«,DY 



+ $a,2$b,2jra 3 



(5.55) 



where Qi is the coefficient of the first order expansion of Q given in (5.E), with /i 2 = in this case and w — > w a 
J 2 m + (D/K eq ) 2 + T 2 /4: 



Qi 



(l + Q^D)m 



J 2 m (3w oo +T/2) 
2«£K +T/2) 



3P„ 



(5.56) 
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and 



a x = 1 + 3£>Q(1 - 2P) + (£>Q) 2 (3 - 4P) + (L>Q) 3 , 
a 2 = 4(1 + £>Q) \kLQ(1 + P) - DQ{1 - P)] , 

a 3 = (1 + DQ - m ) (l + XLq) (5 - 12r + 12r 2 ) . 



(5.57) 
(5.58) 

(5.59) 



In the asymptotic limit this solutions do not coincide with ( 5.39|)- ( 5.42 ). That means that they are also di fferen t 
from the static limit of the correlation functions found in section |2j from the inve rse o f the Hessian matrix (2.26). 
This is due to the fact that the static does not take into account the constraint ( |2.9| ) on the configuration space. 
Above the Kauzmann temperature the dynamics never reaches this constraint so that, even if it is slowed down by 
its existence, it brings to the same static results. But as soon as we perform the dynamics at To or below it, the 
asymptotic regime will never coincide with the equilibrium one. The system will be stucked forever in one ergodic 
component of the phase space, artificially created by imposing the constraint in the dynamics of the model, but not in 
the Hamiltonian. The implementation of the constraint in the dynamics ma kes the variance A 2 of the distribution of 
the random variables giving the updating of Monte Carlo dynamics (see (3.6)) diverging, when T < T Q . The divergent 
factor of A 2 is the quantity T(t) appearing in the equations of motion so far discussed. Going from a regime where 
the contribution of 0(1) are relevant (T > Tq) to another one where they are not even sub-leading, with respect to 
O(r), we loose the static limit. 



We find the solutions of the equations (5.25) for the two-time correlation functions following exactly the approach 
shown in the preceding section, with the following expressions for the functions mi and /: 



mi =4I>iT + 0(/4T) 

/ = -4TT(1 + QD) - 8T -^-— (1 - 3r + 2r 2 ) + 0(^ 2 T) 



1 + QD 

The two-time correlation functions come out to be: 



C lb (t,t') 



r\]b-l r 



[2mi(Q] 
1 + 



x{^J t [l + Q(t")D]T(t")T(t")dt"Y 6=1,2. 
C 2b (t,t') ~2{m 1 (t)] b - 1 C lb (t,t') 



1,2. 



For what concerns the response functions, from (5.33)-( 5.35) we get 



G 12 (t,t + ) 



G 22 (t,i+) 



4Tr 2T(1 - 2r) 



K 



K 



8T(r Ml r 



8miTr 4miT(l-2r) 2 16miT (I>i) 2 16rI>iT 



K 



K 



K 



16m 2 Tr 8m 2 T(l-2r) 2 32m 2 T(I>i) 2 64mirl>iT 8m (l - 2r) 2 T 



K 



K 



K 



K 



where this time the contributions T/ii and (r/xi) 2 are both of order T and we take them into account. 
The two-time behaviour of the response functions is as in ( p5| ) with / given by ( |5.61| ). 
The last thing that we need, before computing the T FD , is the derivative 



d t ,C 11 (t,t') = ^d v Cn(t',t')-fc>)^C n (t,t') 
h(t) h[t) 



4T(i') 



1 



1 + Qoo-D 



rrioQxD 



[1 + Q(t')D]T(t') x 



1 + Qoo-D^ 

exp LtJ [1 + Q(t)D] T(t")T(t") dt' 



(5.60) 
(5.61) 

(5.62) 

(5.63) 

(5.64) 
(5.65) 
(5.66) 



(5.67) 



It follows 
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-iFD (± J 



(t,t')~T e (t') 



1 + 



T e FD (t') 



(5.68) 



In this case 0(1/T) = O(pi) is always smaller than 0(112) because 7 > 1: in the long time regime Tf D (t) coincides 
with T e (t). 



3. Effective temperature from the fluctuation formula 



A self consistent picture with an effective temperature should also imply that the same effective temperature also 
governs other physical variables. From the expression of rrii(t;T) as function of H we can compute the quantity 
^(/O = §m± 1 ihox is the contribution to susceptibility in a cooling-heating setup to a change in the field H at fixed 

time (also called fluctuation susceptibility) . In a cooling experiment the whole susceptibility can, indeed, be written 
as 1 



Xab 



dm a 



dH h 



dm a 



dH h 



+ 



dm a 



T,T C 



dT K 



dT R 



T,H h 



dH h 



dm a 




dm a 


dT e 




dm a 


dT e 




T,t 


dT e 


T,H b dH b 


H 

T,t 


dT e 


T.H b dHb 



x s l uct it) 



xT(t) 



x a T f (t) 



(5.69) 

(5.70) 
(5.71) 



Here we are considering an aging situation, so only the first term is relevant. We can reasonably assume that the 

flu 
A Q 6 



vct (t) can take the form 



xl l b uct (t) 



dm a 



dH h 



= N 



(5m a (i)5m b (t)) fa 



fast 



T.l 



T 



N 



(5m a (t)Sm b (t)) 



T [ J l \t) 



(5.72) 



where (. . ■) f ast / s i ow is the average, respectively, over fast and slow processes. The fast ones being governed by the 
heat-bath temperature, the slow ones by some effective temperature tJ depending on the time scale t. Through 
xi l b UCt (t) one can l°°k a t the connection between the fluctuation effective temperature T e , introduced in [B, and 
the other effective temperatures so far defined. To work it out we start from: 



x{[ uct {t) 



dH 



N 



(<5mi(t)<5mi(f)) Cn(t,t) 



T.l 



T, 



ifi) 



T, 



(fi) 



Using the following expression of toi got from (3.13): 

mi(t;T, H) 



L 
1 



D 



JK (mi(i; T, H),m 2 {t; T, H);T) 



the fluctuation susceptibility xi 



if l) 



turns out to be: 
dm 1 



dH 



T .t K(1 + QD) 



0( Pl ). 



(5.73) 



(5.74) 



(5.75) 



Here we are neglecting terms like dni/dH and d^/dH, of or der h i or hi gher (we deal with the regimes [T > Tq,Vj] 
and [T < To, 7 > 1] where [i\ \i2)- Taking the expressions ( 5.22| ) and ( ^.55| ) we see that in both dynamic regimes 
that we are considering the leading term of C\\ can be written as 



C n (t,t) = 



1 + QD 



o{nl) 



and this leads to 



T, 



(fi) - 



K(m +fJ, 2 ) + O(i^), 



(5.76) 



(5.77) 



thus coinciding with (4.3) to the order of our interest, i. e. 0(^2)- At higher orders there will be non-universalities. 
If 7 < 1 the t erms of 0(^2) become dominant with respect to 0(fi2), leading to the same situation that we had for 
Tf D in ( 5.38| ), namely the thermodynamic description does not lead to a unique effective parameter. 
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6. CONCLUSIONS 



In this paper we consider a model that owns all the basic properties of a fragile glass, built by processes evolving 
on two well separated time scales, representing the a and (3 processes taking place in real glassy materials. 

Also, the model is provided with a constraint applied to the harmonic oscillators dynamics, i. e. to the slow 
processes dynamics, in order to reproduce the behaviour of a good glass former. 

Introducing a particular Monte Carlo dynamics Jl^,|2^,§| and developing it analytically, thus having the opportunity 
of probing it in more detail with respect to a numeric study, we found equations of motion that are in all respect 
those typical of glass relaxation. 

By means of the constrained dynamics we defined the Kauzmann temperature To as the one at which the constraint 
is reached, asymptotically, for the first time in a cooling experiment from high temperature. There we showed how 
the real thermodynamic phase transition JTT[ , taking place due to the breaking of the ergodicity in the landscape of 
our model, rich of degenerate minima, is characterized. 

A detailed study of the dynamics was performed both above and below the Kauzmann temperature and for 
arbitrary values of the exponent 7 generalizing the typical VFTH behaviour, usually assumed for glasses, to 
r eq = exp [A/(T — T )] 7 . The dynamics in the aging regime of both one time and two-time variables has been 
carefully analyzed, including the corrections to this regime, relevant at shorter times. 

The decoupling of time scales is foundamental for a generalization of equilibrium thermodynamics to systems far 
from equilibrium. We tested on our exactly solvable model whether or not the generalized approach holds, involving 
one extra variable, namely the effective temperature, in the description of the non-equilibrium thermodynamics. By 
effective temperature we mean a thermodynamic quantity that would be the temperature of a system at equilibrium 
visiting with the same frequency the same states that the real - out of equilibrium - system at temperature T is 
visiting on a given time-scale during its dynamics. This kind of parameter appears in the thermodynamic functions 
together with the heat-bath temperature and the fields coupled to the system's observables and is coupled to the 
configurational entropy. In our work it has been first derived as the function of time (for given values of the heat 
bath temperature and of all the other parameters of the model) such that the evolving system out of equilibrium can 
be characterized by a probability measure of the configurations having a Boltzmann-Gibbs form with a factor 1/T e 
instead of 1/T in front of the Hamiltonian. 

Generally speaking, in order to recast the out of equilibrium dynamics, the history of a system that is far from 
equilibrium can be expressed by more than one effective parameter. This happens when more than one time-scale are 
involved in the dynamic evolution of a system. In those cases to every time-sector will correspond an effective tem- 
perature |m| . Moreover the number of effective parameters needed to make such a translation into a thermodynamic 
view can, in principle, be equal to the number of relevant observables considered in every time sector. For certain 
dynamic regimes, determined by the temperature and by the VFTH exponent 7, however, the effective parameters 
pertaining to processes having the same time scale become equal to each other in time, for large times. 

As we saw in section ^, in the dynamic regimes reported in section |^, when the variation /ii(t) with respect to 
equilibrium of the variable mi is much greater than the distance from the constraint, ^{t), the effective temperature 
alone is enough for a complete thermodynamic description of the dominant physic phenomena (the effective field 
H e = H) , while in the regimes where [i\ (i) is no more negligible with respect to /12 (t) , the effective field H e (t) is also 
needed to map the dynamics on long time scales into a thermodynamic frame. 

From the time behaviour of the slow varying observables in the aging regime we found in section ^ a VFTH 
relaxation time dependence on temperature above the Kauzmann transition and we derived the Adam-Gibbs relation 
between the relaxation time and the configurational entropy, that we can explicitly compute for our model. 

We have been also able to study the dynamics of the system quenched to a temperature below the Kauzmann 
temperature. At long, finite time t we see that it is possible to introduce an instantaneous relaxation time depending 
on the heat bath temperature in a non-trivial way but expressible in terms of the effective temperature. What we 
got in this way is actually a VFTH law where the heat bath temperature has been substituted by a time dependent 
effective temperature T e (t) and its asymptotic value T e takes the place of the Kauzmann temperature To. Such a 
relation for the time scale of the aging dynamics below To could hold very well in more general systems. 

At equilibrium the heat bath temperature enters many relations that can be rigorously proved and connected to each 
other in the frame of thermodynamics. Out of equilibrium we miss first principles to start with in the generalization of 
such a construction. We do not have any guaranty, for instance, that a given definition of effective temperature, done 
generalizing a given equilibrium formula, would match any other definition coming from the generalization of another 
equilibrium formula. At equilibrium the heat bath temperature enters the Boltzmann-Gibbs measure, the laws of 
thermodynamics, the fluctuation-dissipation theorem, different Maxwell relations, etc. However, out of equilibrium 
we have to check whether a single definition of effective temperature is compatible with any other. Since all effective 
temperatures have definite limits for long times, we should verify that these limits are identical (which happens always) 
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and that the leading approach to these limits coincide (which happens only for 7 > 1). If that works out we find a way 
to completely recast the long time domain of the out of equilibrium dynamics into the language of thermodynamics, 
on a given time sector, of course, well separated from the other time sectors of the glassy dynamics. This behaviour 
may occur if the aging is so slow that the system has time enough to clearly set out an effective temperature before 
going to a lower value of it. 

With this aim we also rederived the effective temperature from the fluctuation dissipation ratio and from the 
fluctuation formula connecting the susceptibility with the fluctuations of the slow variables of the system, namely 
Xj. In section ^]we showed that the effective temperature Tf D defined as the fluctuation dissipation ratio tends 
to the effective temperature T e that we got by the quasi static approach only if ^ is negligible with respect to fi2- 
However, this is true, if 7 is greater than one. Otherwise the corrections of order \T 2 are no more sub-leading and 
change significantly the time evolution of T^ D . Already for 7=1, Tf D — > T e only in the infinite time limit, i. e. 
for time scales longer than those of the considered aging regime. Even above the Kauzmann temperature the value 
of the exponent 7 discriminates between different regimes. For 7 > 1 an out of equilibrium thermodynamics can be 
built in terms of the single additional effective parameter T e . For 7 < 1, T e alone does not give consistent results in 
the generalization of the equilibrium properties to the non equilibrium case. In those cases one also needs an effective 
field H e . However, no universal behaviour for the T e , H e combination has been found. 

We have seen that in case 7 > 1 for temperatures above the Kauzmann temperature both the st atics and the 



dynamics in the aging regime can be described by Vogcl-Fulchcr-Tammann-Hesse laws, see eqs. (4.20), (4.22). Notice 
that these laws are not identical, the static one diverging more strongly, due to a larger prefactor of the divergent term 
in the exponent. We also noticed that in this situation 7 > 1 the aging (well) below the Kauzmann temperature can 
be described in a form very analogous to the aging above it. For this reason it is meaningful to compare experiments 
in those two aging regimes, and in particular, to test whether the decay of measurable quantities, like the energy, the 
volume or the magnetization, has a common temporal law in the full aging regime. 

Finally we checked the consistency of the widely mentioned thermodynamic picture by writing the first and second 
law using the effective temperature and we verified that it can be also computed from the generalization of the Maxwell 
relation at equilibrium giving the heat bath temperature as derivative of internal energy with respect to the entropy 



(see (4.12)), finding the same results we got from the other derivations in the same validity limits. However, we also 
found that, even within the solvable model of this paper, the applicability of the thermodynamic picture may or may 
not be valid, depending on the question how fast the relaxation time diverges. 

ACKNOWLEDGMENTS 

We thank F. Ritort for useful discussions. The research of L. Leuzzi is supported by FOM (The Netherlands). 

APPENDIX: MONTE CARLO INTEGRALS 

Here we present the expressions of the integrals that we use in computing the dynamics of the observables following 
the Monte Carlo method explained in section |j. We recall that x, defined in 3.2, is the energy difference b ewte en 



the current configuration of the system and the one proposed for the updating. The variable r (defined in ( 3.20 ) ) 
is instead the distance of the effective temperature T e from the heat-bath temperature (that is also the equilibrium 
value of T e in the dynamic regime above the Kauzmann temperature). First we define the abbreviation: 



<t:T 

that is the leading term of the acceptance ratio of the Monte Carlo dynamics given by: 



(A.l) 



j dxW(f3x)p(x\m l , ma) = T - - — 2 I±^ + + JL (1 _ 4r + i 6r 2 - 24r 3 + 16r 4 ) + O (J§j \ . (A.2) 



Then we give the behaviour of the derivative with respect to time of the energy 

,.2\ 



J dxW(/3x)xp(x\mi ,m 2 ) = -4rT e T (\ - ^ — 2r + 2r ) + ^_^_ Ur + 2gr 2 _ 3^3 + lQr ^ + Q / J_ 



(A.3) 



2G 



and of the variable mi (defined in ( |2.3| )) 



dxW(f3x)y 1 (x)p(x\m 1 ,m 2 ) = [ T - (1 - 3r + 4r 2 ) + O [ - 



(A.4) 



In section || we compute the correlation and the response functions. In order to find their time dependence we need 
the following derivatives. In these formulae we show the derivatives with respect to mi and m.2, taken as indipendent 
variables, of the effective temperature T e , the variable r and the leading term of the Monte Carlo acceptance ratio T. 
They are: 



dm i 
dT e 
dm? 



2K - mi 

K(P+1): 



(A.5) 
(A.6) 



dr 1 - 3r + 2r 2 
m + H2 
1 - 3r + 2r 2 



9m i 
9m2 



P- - mi ) 



m + ^2 



(P + l); 



and 



9T 

dm i 

or 

dm2 



= -T 
= T 



mi 7 . „ , 1 — 2r 
' (2r+l) + 2 



/'2 



TO + ^2 \ «/ 



I 

P mi 



2r + i (p + i 

2^2 mo + //2 



(A.7) 
(A.8) 

(A.9) 
(A.10) 



Furthermore, we show the extensive computation of the coefficients of equations ( 5.19| )-(5.21) and J5.25 ) for the 
dynamics of the two-time observables: 



fo = 9 mi / dxW{f3x)p(x\mi,m 2 ) ~ -mry- 



M2 



2r - 1 + 4r - 8r 2 + — (13 - 56r + 136r 2 - 160r 3 + 80r 4 ) 



O 



T 



to + fi2 



L 

P 7 -mi 



1 - 2r ■ 
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f 



3 1 



T \ 
fa J 
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2^2 
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92 = d m2 J dxW((3x)y 2 (x)p(x\mi,m 2 ) ~ 2mi#i + -~;g 3 + 8Tr 1 + q D 



QD 
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93 = d m2 / dxW((3x)xp(x\m 1 ,m 2 ) ~-2jTT e 



2r-5+12r-12r 2 --(3 + 9r-54r 2 + 140r 3 -160r 4 +80r 5 ) 



O 



1 - 3r + 2r 2 - - (l - 5r + 12r 2 - 14r 3 + 8r 4 



01 s 



(A.18) 



(A.19) 



All partial derivatives with respect to mi have been computed keeping mi fixed a nd vic e- versa . At this stage time 
has not yet been introduced. Introducing it we are able to make an expansion of (|A~ll)-( |A~T8| ) in powers of ^(t). 
In the formula we already performed such an expansion, breaking it at O (T/(r^2)) that is more than sufficiently 
refined to derive the dynamics of the correlation and response functions in all the regimes of our interest. 
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